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Abstract 

Group based moving frames have a wide range of applications, from the classical equiva- 
lence problems in differential geometry to more modern applications such as computer vision. 
Here we describe what we call a discrete group based moving frame, which is essentially a 
sequence of moving frames with overlapping domains. We demonstrate a small set of gen- 
erators of the algebra of invariants, which we call the discrete Maurer-Cartan invariants, 
for which there are recursion formulae. We show that this offers significant computational 
advantages over a single moving frame for our study of discrete integrable systems. We 
demonstrate that the discrete analogues of some curvature flows lead naturally to Hamilto- 
nian pairs, which generate integrable differential-difference systems. In particular, we show 
that in the centro-afhne plane and the projective space, the Hamiltonian pairs obtained can 
be transformed into the known Hamiltonian pairs for the Toda and modified Volterra lattices 
respectively under Miura transformations. We also show that a specified invariant map of 
polygons in the centro-affine plane can be transformed to the integrable discretization of the 
Toda Lattice. Moreover, we describe in detail the case of discrete flows in the homogeneous 
2-sphere and we obtain realizations of equations of Volterra type as evolutions of polygons 
on the sphere. 

Dedicated to Peter Olver in celebration of his 60th birthday 

1 Introduction 

The notion of a moving frame is associated with Elie Cartan [3j, who used it to solve equiv- 
alence problems in differential geometry. Moving frames were further developed and applied 
in a substantial body of work, in particular to differential geometry and (exterior) differential 
systems, see for example papers by Green and Griffiths |12j . From the point of view of 
symbolic computation, a breakthrough in the understanding of Cartan's methods came in a 
series of papers by Fels and Olver [61 [7J , Olver [37] , Hubert [151 EE E] , and Hubert and 
Kogan \18\ I19j. which provide a coherent, rigorous and constructive moving frame method free 
from any particular application, and hence applicable to a huge range of examples, from classical 
invariant theory to numerical schemes. 

For the study of differential invariants, one of the main results of the Fels and Olver papers is the 
derivation of symbolic formulae for differential invariants and their invariant differentiation. The 
book |29j contains a detailed exposition of the calculations for the resulting symbolic invariant 
calculus. Applications include the integration of Lie group invariant differential equations, to 
the Calculus of Variations and Noether's Theorem, (see also |25t HOj). and to integrable systems 

( [23 GEl ESSES). 
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The first results for the computation of discrete invariants using group-based moving frames 
were given by Olver [36J who calls them joint invariants; modern applications to date include 
computer vision [37] and numerical schemes for systems with a Lie symmetry [3J 122] [23], 123] [30] . 
While moving frames for discrete applications as formulated by Olver do give generating sets of 
discrete invariants, the recursion formulae for differential invariants that were so successful for 
the application of moving frames to calculus based applications do not generalize well to these 
discrete invariants. In particular, these generators do not seem to have recursion formulae under 
the shift operator that are computationally useful. To overcome this computational problem, 
we introduce a discrete moving frame which is essentially a sequence of frame^j and prove 
discrete recursion formulae for a small computable generating sets of invariants, which we call 
the discrete Maurer-Cartan invariants. 

We show that our definitions and constructions arise naturally and are useful for the study of 
discrete integrable systems. These arise as analogues of curvature flows for polygon evolutions in 
homogeneous spaces, and this is the focus of the second half of the paper. The study of discrete 
integrable systems is rather new. It began with discretising continuous integrable systems in 
1970s. The most well known discretization of the Korteweg-de Vries equation (KdV) is the Toda 
lattice [30] 

d 2 u s 

= exp(n s _i - u s ) - exp(u s - u s+1 ). (1) 

Here the dependent variable u is a function of time t and discrete variable s £ Z. We can obtain 
a finite-dimensional version by picking N € N and restricting to 1 < s < N subject to one of 
two types of boundary conditions: open-end («o = un = 0) or periodic (u s+ n = u s for all s and 
some period N). Using the Flaschka [8, 9J coordinates 

q s = -jp Ps = expf/u, - u s+ i), 

we rewrite the Toda lattice in the form 

dps , s dq s 

-^r = Ps{q s ~ q.s+i), — =p s -i-p s . (2) 

Its complete integrability was first established by Flaschka and Manakov [8J [9j [28] . They con- 
structed the Lax representation of system ^ and further solved it by the inverse scattering 
method. 

Another famous integrable discretization of the KdV equation is the Volterra lattice [28[ [2"T] 

dq s , . 
— = q s {q s+ i - q s -i). 

By the Miura transformation q s = p s Ps—i, it is related to the equation 

P 2 s (Ps+i -Ps-i), (3) 



dPs 2, 



dt 

which is the modified Volterra lattice, an integrable discretization of the modified KdV equation. 

Since the establishment of their integrability, a great deal of work has been contributed to the 
study of their other integrable properties including Hamiltonian structures, higher symmetry 



1 A sequence of moving frames was also used in [23] to minimize the accumulation of errors in an invariant 
numerical method. 
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flows and r-matrix structures, as well as to the establishment of integrability for other systems 
and further discretising differential-difference integrable systems to obtain integrable maps. For 
example, the time discretisation of the Toda lattice ^ leads to the integrable map (p, q) h-> (p, q) 
defined by 

Ps+i ~ . f Ps Ps-l\ , A s 

Ps=Ps ls7> qs = qs + c \J s -J^)' (4) 

where c € M is constant and the function f3 is given by the recurrent relation 

p s = 1 + cq s - c 



Ps-l 

Some historical background about the development of the theory of discrete integrable systems 
can be found in [39] ■ Some classification results for such integrable systems including the Toda 
and Volterra lattice were obtained by the symmetry approach 



In this paper, we introduce the concept of discrete moving frames, and under conditions which are 
satisfied for the range of examples we study, we prove theorems analogous to the classical results 
of the continuous case: generating properties of Maurer-Cartan invariants, a replacement rule, 
recursion formulas, and general formulas for invariant evolutions of polygons. Once the ground 
work is in place, we study the evolution induced on the Maurer-Cartan invariants by invariant 
evolutions of iV-gons, the so-called invariantizations. We consider the resulting equations to 
be defined on infinite lattices, i.e. s £ 7L for both iV-periodic and non-periodic cases. We will 
show that the invariantization of certain time evolutions of iV-gons (or so-called twisted N-gons 
in the periodic case) in the centro-affine plane and the projective line MP 1 naturally lead to 
Hamiltonian pairs. Under the Miura transformations, we can transform the Hamiltonian pairs 
into the known Hamiltonian pairs for the Toda lattice (|2j) pQ and modified Volterra lattice Q 
|26j respectively. We also show that a specified invariant map of polygons in the centro-affine 
plane naturally leads to the integrable map Q for the time discretisation of the Toda Lattice. 
We will analyze in detail the case of the 2-homogeneous sphere. We will use normalization 
equations to obtain Maurer-Cartan invariants and we will prove that they are the classical 
discrete arc-lengths (the length of the arcs joining vertices) and the discrete curvatures (tt minus 
the angle between two consecutive sides of the polygon) . We will then write the general formula 
for invariant evolutions of polygons on the sphere and their invariantizations. We finally find 
an evolution of polygons whose invariantization is a completely integrable evolution of Volterra 
type. 

The arrangement of the paper is as follows: In Section 2 we introduce moving frames and the 
definitions and calculations we will discretize. In Section 3 we introduce discrete moving frames 
and discrete Maurer-Cartan invariants and prove the main Theorems of the first part of the pa- 
per concerning these invariants. In Section 4 we begin our study of Lie group invariant discrete 
evolutions of iV-gons in a homogeneous space. We show the discrete moving frame yields an 
effective and straightforward reduction or invariantization of the evolution, to produce discrete 
analogues of curvature flows. In Section 5 we describe discrete invariant evolutions and demon- 
strate their integrability. In Sections 5.1 and 5.2 we study the centro-affine and the projective 
cases, including the completely integrable evolutions of polygons and their associated biHamilto- 
nian pair. We also briefly describe invariant maps and their invariantizations and demonstrate 
that the invariant map in the centro-affine plane leads to the integrable discretization of the 
Toda lattice ([2]). However, we leave their thorough study, including a geometric interpretation 
of integrable maps and more examples of maps with a biPoisson invariantization, for a later 
paper. Finally, Sections 5.3 describes the more involved example, that of discrete evolutions on 
the homogeneous sphere. We conclude with indications of future work. 



3 



2 Background and definitions 



We assume a smooth Lie group action on a manifold M given byGxM-f M. In the significant 
examples we discuss in the later sections, M will be the set of iV-gons in a homogeneous space, 
that is, iV-gons in G/H where H is a closed Lie subgroup, or more generally M = (G/H) with 
the standard action. But in what follows this restriction does not need to be in place. 

2.1 Moving Frames 

We begin with an action of a Lie Group G on a manifold M. 

Definition 2.1. A group action of G on M is a map G x M — > M, written as (g, z) i— > g ■ z, 
which satisfies either g ■ (h- z) = (gh) ■ z, called a left action, or g ■ (h- z) = (hg) ■ z, called a right 
action. 

Definition 2.2 (Invariants). Given a smooth Lie group action G x M — > M, a function / : 
M — > K is an invariant of the action if I(g ■ z) = I(z) for any g £ G and any z € M. 

The set of invariants form an algebra; here we consider local invariants and they will typically 
be locally smooth. When we talk about generators of the algebra we are referring to functional 
generators. 

Let us write g ■ z as z to ease the exposition in places. Further, we assume the action is free 
and regular in some domain Q C M, which means, in effect, that for every x £ Q there is a 
neighbourhood U C of x such that: 

1. the intersection of the orbits with U have the dimension of the group G and further foliate 
U; 

2. there exists a submanifold /C C U that intersects the orbits of U transversally, and the 
intersection of an orbit of IA with fC is a single point. This submanifold K, is known as the 
cross-section and has dimension equal to dim(M) — dim(G); 

3. if we let O(z) denote the orbit through z, then the element h £ G that takes z E U to k, 
where {k} = 0(z) n /C, is unique. 

Under these conditions, we can make the following definitions. 

Definition 2.3 (Moving frame). Given a smooth Lie group action G x M — > M, a moving 
frame is an equivariant map p : U C M — > G. We say U is the domain of the frame. 

Given a cross-section fC to the orbits of a free and regular action, we can define the map p :U — >• G 
such that p(z) is the unique element in G which satisfies 

p(z)-z = k, {k} = O(z) n K, 

see Figure [TJ We say p is the right moving frame relative to the cross-section /C, and fC provides 
the normalization of p. This process is familiar to many readers: it is well known that if we 
translate a planar curve so that a point p in the curve is moved to the origin, and we rotate 
it so that the curve is tangent to the x-axis, the second term in the Taylor expansion at p 
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orbits 



Figure 1: The definition of a right moving frame for a free and regular group action. It can be 
seen that p(g • z) = p{z)g~ l (for a left action). A left moving frame is obtained by taking the 
inverse of p(z). 

is the Euclidean curvature at p. The element of the Euclidean group taking the curve to its 
normalization is indeed a right moving frame. 

By construction, we have for a left action and a right moving frame that p(g • z) = p{z)g~ l so 
that p is indeed equivariant. A left moving frame is the inverse of a right moving frame, so for 
a left action and a left moving frame, the equivariance is p{g ■ z) = gp{z). The cross-section 
K, is not unique, and is usually selected to simplify the calculations for a given application. 
Typically, moving frames exist only locally, that is, in some open domain in M. In what follows, 
we conflate IA with M to ease the exposition. In applications however, the choice of domain may 
be critical. 

In practice, the procedure to find a right moving frame is as follows: 

1. define the cross-section K, to be the locus of the set of equations tpi(z) = 0, for i = 1, r, 
where r is the dimension of the group G; 

2. find the group element in G which maps z to k G fC by solving the normalization equations, 



Hence, the frame p satisfies ipi(p(z) ■ z) = 0, i = 1, r. 

Remark 2.4. In practice, both left and right actions occur naturally, as do left and right moving 
frames, and the calculations can be considerably easier in one handedness than the other, so 
both occur in the examples. However, since the handedness can be changed by taking inverses, 
there is no need to consider both in a theoretical development. We will usually work with right 
and left moving frames for a left action. 

In what follows we will use the terms right and left frames to designate left and right moving 
frames since this creates no confusion. Invariants of the group action are easily obtained. 

Definition 2.5 (Normalized invariants). Given a left or right action G x M — > M and a right 
frame p, the normalized invariants are the coordinates of I(z) = p{z) ■ z. 

Indeed, for a left action we have 



i/ji(z) = ipi(g • z) = 0, 



i = 1, ...,r. 



I{9 ■ z) = p{g ■ z) ■ g ■ z = p(z)g g-z = p(z) ■ z 



I{z). 
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The normalized invariants are important because any invariant can be written in terms of them; 
this follows from the following: 



Theorem 2.6 (Replacement Rule). If F(z) is an invariant of the action G x M — > M ; and 
I(z) is the normalized invariant for a moving frame p on M, then F(z) = F(I(z)). 

The theorem is proved by noting that (for a right moving frame) F(z) = F(g ■ z) = F(p(z) ■ z) = 
F(I{z)) where the first equality holds for all g G G as F is invariant, the second by virtue of 
setting g = p(z), and the third by the definition of I(z). 

The Replacement Rule shows that the normalized invariants form a set of generators for the 
algebra of invariants. Further, if we know a sufficient number of invariants, for example, they 
may be known historically or through physical considerations, then the Replacement Rule allows 
us to calculate the normalized invariants without knowing the frame. 

We can apply this theory to product (also called diagonal) actions, the action with which our 
paper is concerned. This move has its shortcomings, as we see next. 

Given a Lie group action G x M — > M, (g, z) i— )■ g ■ z, the product action is 



G x (M x M x • • • x M), (g, (z 1 ,z 2 , . . -,z N )) (->• (g ■ zi, g ■ z 2 , ■ ■ ■ , g ■ z N ). 
In this case M = M N = M x • • • x M and the normalized invariants are the invariantized 



These were called joint invariants in [36J . 

Remark 2.7. From now on the manifold where G acts will be the product M = M x M ■ ■ ■ x 
M, hence, questions like freedom of the action, etc, will refer to the diagonal action on the 
product. We will also assume that the number of copies N of the manifold M is high enough to 
guarantee that the action is free. In fact, there is a minimal number that achieves this (called 
the stabilization order) , the proof of this result can be found in [2] . 

The following simple scaling and translation group action on K will be developed as our main 
expository example, as the calculations are easily seen. 

Example 2.8. 

Let G = R + x R = {(A, a) |, A > 0, a e K} act on M = R as 



The product action is then given by z n h-» Xz n + a for all n = 1, . . . , N. We write this action as 
a left linear action as 



2 In the literature the invariantization of a component z t is sometimes denoted as I(zi), i(zi) or t(zj), here we 
use Ii(z) 




F(z 1 ,z 2 , . . . ,z N )= F(I 1 (z),I 2 (z), . . .,I N (z)) 



z i— )■ Xz + a. 
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There are two group parameters and so we need two independent normalization equations. We 
may set 

9 ■ z\ = 0, g • z 2 = 1 

which are solvable if z\ 7^ Z2, and this then defines the domain of this frame. In matrix form 
the right frame is, 

1 zi 

P(Z1,Z2, ■ ■ ■ , Zn) = I Zl — 22 Z1 — Z2 

1 



Equivariance is easily shown. For h = h(n,b): 



1 [izi + 6 



p(h ■ zi, h ■ Z2, ■ ■ ■ , h ■ zn) = \ fi(zi - z 2 ) fi(zi - z 2 ) 

1 

1 zi 



Zl - Z 2 Zl- Z 2 

1 

= p(zi,Z2,...,z N )h(n,b)' 
The normalized invariants are then the components of 

1 Zl 









II 


fi 


i) 




V 




1 







zi - Zl Zl- z 2 
1 



Zl z 2 z 3 ■■■ z N 

1 1 1 ••• 1 



h(z) I 2 (z) h(z) ••• I N (z) 

1 1 1 ••• 1 

Z3 - Zl Z N - Zl 



1 

z 2 - Zl Zl - Zl 

1 1 1 1 

noting that h(z) = and I 2 (z) = 1 are the normalization equations. Any other invariant can 
be written in terms of these by the Replacement Rule. For example, it is easily verified that the 
invariant 



Z5 ~ Z4 \ ( h(z) - h(z) 



zn-zqJ \In{z) - h{z 



Z5 ~ 


Zl 


Z4 ~ 


Zl 


\ 


z 2 - 


Zl 


z 2 - 


Zl 




ZN ~ 


Zl 


Z6 ~ 


Zl 




z 2 - 


Zl 


z 2 - 


Zl 


) 



The above example shows both the power and the limits of the moving frame. In applica- 
tions involving discrete systems, use of the shift operator T taking z n to z n +i is central to 
the calculations and formulae involved. However, it can be seen in the above example that 
T{I n {z)) 7^ I n +i(z) and that instead the Replacement Rule will lead to expressions for T(I n (z)) 
that are as complicated in their new formulae as in the original co-ordinates. By contrast, our 
discrete moving frame theory yields a single expression which, together with its shifts, generates 
all invariants. 

As motivation for our discrete moving frame, the next example will show how an indexed sequence 
of moving frames arises naturally in a discrete variational problem. This sequence of frames 
arises when Noether's conservation laws for a discrete variational problem with a Lie group 
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symmetry, is written, as far as possible, in terms of invariants. As for the smooth case, the non- 
invariant part turns out to be equivariant with respect to the group action, and so by definition 
forms a moving frame, indexed by n in this discrete case. These frames are not constructed 
by normalisation equations and the calculation of the conservation laws makes no use of the 
concepts of moving frame theory. Nevertheless, it is fascinating that frames appear, calculated 
"for free" by the formulae for the laws. 

Example 2.9. 

Consider the discrete variational problem 

for which it is desired to find the sequence (z n ) which minimises C[z], possibly subject to 
certain boundary conditions. We note that J n+ k = ( z n+k+2 — z n+k+i) I '(^n+fc+i — -Zn+fc) is 



invariant under the scaling and translation action, z n i— > Xz n + a (see Example 2.8) and that 
T(J n +k) = Jn+k+i where T is the shift operator. The discrete Euler Lagrange equation for a 
second order Lagrangian in a single discrete variable is (cf. |20| ) 







dL r , 



+ T 



.1 dL r 



+ T- 



dL r 



dz n dz n+ i dz n+2 

which for this example can, after some calculation, be written as 

n — T 2 — T 2 — _l t T 2 

u — J n+1 ,J n ,J n ' °n d n—l 

after shifting to balance the indices. The Lie group invariance of the summand L n implies there 
are two conservation laws (in this case, first integrals) arising from the discrete analogue of 
Noether's Theorem |14} 120] . The general formulae for these are complicated and we do not 
record them here. The two first integrals can be written in matrix form as 



ci 

C2 



^ n (z)v n (J) 



/ 



V 



1 



Zn+1 Z n 
Z n +1 







1 



Zn+1 



\ 



( Jn+1 \ 



\ 



J ft 



rJ <n 



J 



where this defines the matrix A n (z) and the vector of invariants v n (J), and the Cj are the 
constants of integration. It can be seen that in this example, the Euler Lagrange equation gives 
a recurrence relation for J n and that once this is solved, the conservation laws yield (if c\ ^ 0), 



z n+l 



The matrix ^4 n (z) is equivariant under the group action, for each n, indeed 



A n (g ■ z) 



/ 1 



\ 



Zn) 

Zn+1 + b/n 
Zn+1 Z n 



\ 



J 





( 


1 




1 
















b 








V 




l J 


\ 



1 



Zn+1 Z r 
Zn+1 



Zn+1 



Z'ii 



and hence each A n (z) is a left moving frame for the group G = M + k M, albeit for a different 
representation for this group. In fact, this is the Adjoint representation of G, placing this result 
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Figure 2: The location of the arguments occurring in the normalization equations for a moving 
frame for the product action on M x M x • • • x M, shown here as disjoint copies of M for clarity, 
for the Example |2.8[ 




Figure 3: The sequence of moving frames for the product action on M x M x • • • x M in 
Example 2.9 shown here as disjoint copies of M so the location of the arguments occurring 
in the normalization equations can be easily seen. For this example, each /C, is a shift of the 
previous. This sequence of moving frames is an example of a discrete moving frame. 



in line with Theorems on the equivariance of Noether's conservation laws for smooth systems 

ma eh]. 

In Figure [2] is shown, schematically, the location of the arguments occurring in the normalization 
equations for the frame of Example 2.8. By contrast, in Figure [3] is shown, schematically, the 



location of the arguments occurring in the normalization equations for the sequence of frames 



calculated in Example 2.9 In fact, for that example, each of the /Q will be shifts of each other. 



In the next section, we will define a discrete moving frame to be a sequence of moving frames 
with a nontrivial intersection of domains, and will explore the properties of the structure of the 
algebra of invariants that arises. 
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3 Discrete moving frames and the example of twisted iV-gons in 
a homogeneous manifold G/H 



3.1 Discrete moving frames and discrete invariants 



In this section we will state the definition of discrete group-based moving frame along A-gons 
and show that a parallel normalization process will produce right discrete moving frames. A 
discrete moving frame gives many sets of generators for the algebra of invariants under the 
action of the group. Under conditions that are satisfied in our examples, we will show recursion 
relations between these and find a small, useful set of generators which we will denote as discrete 
Maurer-Cartan invariants. 

Our next definition represents the discrete analog of the group-based moving frame we described 
in Section 2. 

Definition 3.1 (Discrete moving frame). Let G N denote the Cartesian product of N copies of 
the group G. Allow G to act on the left on G N using the diagonal action g ■ (g r ) = (gg r ). We 
also consider what we call the "right inverse action" g ■ (g r ) = {gr-g -1 )- (We note these are both 



left actions according to Definition (2.1).) We say a map 

p:M N ^G N 

is a left (resp. right) discrete moving frame if p is equivariant with respect to the diagonal 
action of G on M N and the left (resp. right inverse) diagonal action of G on G N . Since 
p((x r )) G G N , we will denote by p s its sth component, that is p = (p s ), where p s ((x r )) G G for 
all s. Equi variance means, 

Ps(g • {x r )) = p s ((g • x r )) = gp s {(x r )) (resp. p s {{x r ))g~ l ) 

for every s. Clearly, if p = (p s ) is a left moving frame, then p = (p^ 1 ) is a right moving frame. 

Definition 3.2. Let F : M N — > M be a function defined on A-gons. We say that F is a discrete 
invariant if 

F((g ■ x r )) = F((x r )) (5) 

for any g G G and any (x r ) G M N . 



Notice that the quantities p s ({x r )) ■ = II are always invariant as we can readily see from 
9 • I k = Ps(g • (av)) • {9 ■ x k ) = p s {{x r ))g~ l g ■ x k = p s {{x r )) ■ x k = I s k . 

Proposition 3.3 (s-Replacement rule). If p is a right moving frame, and F((x r )) is any in- 
variant, then 

F({x r )) = F({p s -x r )) (6) 

Further, the invariants p s ■ x r = I® with s fixed, r = 1, ...JV generate all other discrete 
invariants. We call the I£ with s fixed, the s-basic invariants. 



Proof. This is the same as Theorem 2.6 for each s. From ([5]) we need only choose g = p s in ([5]). 
In particular, Equation ^ means F can be written as a function of the s-basic invariants: one 
merely needs to substitute x r by p s ■ x r for s fixed. For those familiar with moving frames, this 
proposition follows directly from the fact that each p s is a moving frame and p is a collection of 
N moving frames. □ 
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Locally, a discrete moving frame is uniquely determined by the choice of cross-sections /C s to 
the group orbit through (x r ), as in figure [3j Often these sections are shifts of a first section, 
but that need not be the case in general. The right moving frame component p s is the unique 
element of the group that takes (x r ) to the cross section IC S . Notice that p s will depend, in 
general, on several different coordinates, not just x s . 

Theorem 3.4. Let fC be N local choices /Ci, . . . /Cat of cross-sections to the orbit of G through 
(x r ). Let p = (p s ) 6 G be uniquely determined by the condition 

Ps ■ (x r ) G fC s , (7) 

for any s. Then p = (p s ((x r ))) is a right moving frame along the N-gon (x r ). 



Proof. The proof is entirely analogous to that of a single frame. We will denote by p s ((x r )) 
the unique element of the group determined by equations Q. On the other hand p s ((x r ))(7 _1 
also satisfies those equations, and so p s (g • (x r )) = p s ((x r ))g . Therefore p is a right discrete 
moving frame. □ 



As in the continuous case, a practical way to obtain a right discrete moving frame is using 
a coordinate cross-section, defined as follows. Denote by x l the coordinates of x G M. Let 
p = (p s ((x r ))) G G be uniquely determined by a number of equations of the form 

(p s (( Xr )) ■ X k y = (4)\ (8) 

for any s, k in a set that depends on s, and some chosen values of i that depend on s and k, 
and where (cf,) 1 G R. From now on those (x r ) for which the map g — > (g ■ x r ) 1 has full rank for 
the choice of i's defining the moving frame will be called regular points. We usually refer to the 
{c s k ) 1 G M as the normalization constants. 

Remark 3.5. In some of our examples, and indeed for all examples we envisage, the normal- 
isation equations for the p r are shifts of each other. However, it is interesting that this is not 



necessary to produce the recursion relation (13) we derive in Section 3.2 When normalization 



equations are shifts of each other the construction of a discrete moving frame depends only 
on a single function of several variables and the different components of the moving frame are 
generated by evaluating this function on the corresponding vertices. 



3.2 Discrete Maurer— Cartan invariants 

A discrete moving frame on M N provides ./V different sets of generators, since each p s generates 
a complete set. With this abundance of choices one is hard pressed to chose certain distinguished 
invariants, and a small group that will generate all others. Although a correct choice will depend 
on what you would like to use the invariants for, a good choice in our case are what we will 
denote the discrete Maurer-Cartan invariants. They are produced by the discrete equivalent of 
the Serret-Frenet equations. From now on we will assume that G C GL(m, K) and so K r are 
represented by matrices; this is not necessary, but it is convenient. 

The following definition appeared in |38j. 

Definition 3.6 (Twisted N-gon). A twisted N-gon in a manifold M is a map <j) : Z — > M such 
that for some fixed g G G we have <j)(k + N) = g ■ 4>{k) for all k G Z. (Recall • represents the 
action of G on M given by left multiplication on representatives of the class.) The element 
g G G is called the monodromy of the gon. 
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Remark 3.7. The choice of twisted A^-gons is one among several ways to make the shift operator 
well-defined. Because of the need to apply the shift operator freely on the iV-gons, one is forced 
to either consider infinite polygons, or to impose some kind of periodicity condition. Although 
one could merely choose closed gons, there is an important reason to work with twisted Ar- 
gons instead: twisted A^-gons will give rise to A^-periodic invariants; however, a general set of 
periodic invariants will not, in general, be associated to A^-periodic gons, but to twisted ones. 
Nevertheless, for much of what follows we can also consider infinite gons and we will occasionally 
consider the infinite case in the applications. 



The space of twisted A^-gons in M can be identified with the Cartesian product of A^ copies of 
the manifold M, and hence our previous theory applies. 

Definition 3.8. Let (p s ) be a left (resp. right) discrete moving frame evaluated along a twisted 
N-gon, the element of the group 

K s = pJ 1 P*+i ( res P- Ps+iPs 1 ) 

is called the left (resp. right) s-Maurer-Cartan matrix for p. We will call the equation p s +i = 
p s K s the discrete left s-Serret-Frenet equation. 



One can directly check that if K s is a left Maurer-Cartan matrix for the left frame (p s ), then 
K~ is a right one for the right frame p= (p^ 1 ), and vice versa. 

The equivariance of p immediately yields that the K s are invariant under the action of G. We 
show next how the components of the Maurer-Cartan matrices generate the algebra of invariants 
by exhibiting recursion relations between them and the normalized invariants, for our expository 
example. 



Example 2.8 cont. Recall the group is G = M + ix R acting on M as z = (A, a) ■ z = Xz + a, 











r 

\ 




i 





which is represented so that the action is left (multiplication) as 

Xz + a 
1 

If we take the normalization equations of p s to be z s = 1 and z s+ \ = 0, then the right discrete 
frame is 

p s = ^s+l Zg Z s +i Z s 

V 1 

By definition, p s • z r = 1% = —(z r — z a +i)/(z s +i — z s ) and therj^] 

Ks = PsPjli = ( + 2 J Y 2 ) (9) 

as can be verified directly. We do not need, however, to have solved for the frame to obtain this 
result; it can be obtained directly by noting that we have both 



z s +l z s+2 \ _ ( 1 \ / z s+1 z s+2 \ _ ( /; 



p, +1 [ ~r T) = {ll)- T T ) = {! T ) 

3 We calculate the left Maurer-Cartan matrix for the left frame p _1 as the calculations are simpler. 



rs 
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or 



PsP s+ l 



/; 
1 



g+2 
1 



1 

1 1 



verifying the result in Q. Note that for the second calculation, we are using the invariants I, 
symbolically; if we have not solved for the frame, then we may not know what these are. 



Calculations similar to those of Equation (10) give recursion relations between the I*. For 
example, we have both 



z g+k z s+k+l 
1 1 



TS Tg 

1 g+k ± s+k+l 



1 



1 



so that 



PsPs+k 



Ps+k 



L s+k 



z s+k z s+k+l 
1 1 



1 

1 1 







TS TS 

± s+k+l ± s+k+l 



K s K s +i ■ ■ ■ K s +k-i 
However, we also have by direct calculation using our Ki, for example, 

K s K s+l 



1 



TS TS+1 TS 1 1 _ TS+1\ 

1 s+2 I s+3 1 s+2\ 1 1 s+3j 







1 



and thus for k = 2 we have J|. 3 = Ig+ 2 (^ ~ -^+3)- ^ n this wa Y) it is possible to obtain all the ij? 
from the components of the K^'s, that is, the This smaller set of generating invariants we 

will denote as the discrete Maurer-Cartan invariants. 

Example 3.9. 

Our next example is the centro-amne action of SL(2, R) on IR 2 ; that is SL(2,M) acts linearly on 
M?. We may cast the centro-amne plane in the form G/H by identifying IR 2 with SL(2,M)/H, 

where H is the isotropy subgroup of e2 = f 1 j • There is a reason why we are choosing e2 rather 



than e\\ the reader will readily see the connection to Example |3. 10 the projective line. In order 
to find a moving frame we will use the following normalization equations 



Ps ' 2Cs 



e 2 , p s -x s+1 



-s+l 



where is the determinant of the two 2-vectors; specifically, the second component of 

p s ■ x s +i being zero is the normalization equation, while the first is determined by the previous 
normalizations and the fact that p s G SL(2,K). This results in the left moving frame 



-l 



I Xs+l X. 



\x s ,3? s _|_l I 

From here, a complete set of generating invariants will be the components of 



Ps ' Xr 



PsXr 



\X r , X s 



1 



\xs,x s +l\ 



\X r , 3J S _|_1 



with s fixed and r = 1, ... TV; that is, \x s , x r \, |x s+ i, x r \ with s fixed and r = 1, . . . N. To see 
how the r-basic invariants will generate the A;-basic invariants it suffices to recognise the relation 



\Xri Xk\ 



\p s X r: p s Xk\ 



\Xf , X s I 
\x r , X s -\-i I 



[Xfc, x s \ 
I Xk 1 Xgjj-X I 
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\Xsi ^s + l | 



obtained by the Replacement Rule, Proposition 3.3 The left Maurer-Cartan matrix is given by 



where k 2 = jj^fg^ and A; 1 = | x s ,x s +i|. By calculations similar to those of the previous 
example, we have that these invariants, with s = 1,. ..,N generate all other invariants. This 
follows from observing 



p s ■ X r = p s p s+1 p s+ xp s+2 ■ ■ ■ p r -ip r 1 {p r ■ X r ) = K s K s+ i ■ ■ • if r _ie 2 

and similarly if r < s. In fact, using kl +1 we could simplify these generators to the simpler set 
\x s , x s +i\ , \x s ,x s+ 2\, s = 1, . . . , N, as expected. 

Our final expository example before proving our result concerning the discrete Maurer-Cartan 
invariants is the projective action of SL(2,K) on iV-gons in MP 1 . 

Example 3.10. 

Consider local coordinates in MP 1 such that a lift from RP 1 to M 2 is given by x — > (^j • ^ n ^ na ^ 

case MP 1 can be identified with SL(2,M.)/H, where H is the isotropy subgroup of x = 0. The 
action is given by the fractional transformations 

a., b s \ a s x r + b s 



Ps ' X r — 1 - I • X r — i — , G's^s OsCs — -L" 

\ c s Us J C s X r + d s 

In this particular case one can use a geometric description to find a simple moving frame without 
resorting to normalization equations (most simple normalization equations would produce a 
moving frame with no clear geometric meaning). Indeed, we can lift x s to V s <E M. 2 so that 

det(V s+ i, V s ) = 1 for all s. It suffices to define V s = t s ( Xs ) and solve 



t s t s+1 det 



Xs 3?s+l 
1 1 



If N is not even, this equation can be uniquely solved for t s , s = 1, . . . , N, using the twisted 
condition t^ +s = t s for all s. The element 

p s = (V s+1 , V s ) 

is clearly a moving frame since the lift of the projective action to M 2 is the linear action. It also 
satisfies p s ■ o = x s , where o = (0, 1) T is the equivalence class of H in SL(2,IR)/ff, where recall 
H is the isotropy subgroup of x = 0. Given that M? is generated by V s , V s +± for any choice of s, 
we have that 

V s +2 = k s V s+ i - V s 

for all s (the coefficient of V s reflects the fact that det(Vs + i, V s ) = 1 for all s). From here 

fk s 1 

ps+1 = Ps -1 
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and so the Maurer-Cartan matrix is given by 

*- = (-i J)' (12) 

This invariant appeared in |38j . 

Directly from the definition of Maurer-Cartan matrix one can see that their entries, together 
with Ig, are generators for all discrete invariants of iV-gons. 

Proposition 3.11. Let K = (K s ) be the left (resp. right) Maurer-Cartan matrix associated to 
a moving frame (p s )- Then K satisfies the recursion relations 

K s • J r s+1 = I» {resp. K s ■ I s r = I s r +1 ) (13) 

for all s,r. Furthermore, all discrete invariants are generated by K and J|, s = 1, . . . JV. 



Proof. The proof is a direct straightforward calculation. Since K s = p s p s +i, K s will satisfy 

K ■ T s+1 — rT x ■ n ■ T s — T s 

as stated. Since the right Maurer-Cartan matrix for the right frame p s is the inverse of the left 
one for pj 1 , it will satisfy K s ■ I* = I£ + ■ 

To prove the generating property, let us denote by / any discrete invariant. We know that I will 
be a function of the basic invariants I* by the Replacement Rule, so it suffices to show that 1^, 
r 7^ s, are all generated by the entries of K s and J*. And this directly follows from the recursion 
formula. Indeed, we can use the recursion formulas for -fT s _i, . . . , Kq to find I§ , ■ ■ ■ , ij , I® from 
the values I| recursively. If we invert the recursion equation to read 



rs+l 



K7 1 ■ It 



then we can use the equations for K s , K s+ \, . . . K^-i to equally generate J| +1 , I s s +2 , . . . , 1^ 
from the value I*. Therefore, K s and I£ will generate I* for any s, r = 1, . . . , N. □ 

Remark 3.12. If the normalization equations guarantee that p is uniquely determined, the 
rank of the map g — > g ■ will have maximum rank for each fixed s. Thus, it will equally 



guarantee that K = (K s ) is uniquely determined by the recursion relations (13) directly from 
the normalization constants. Therefore, one is able to find the Maurer-Cartan matrix without 
knowing explicitly the moving frame, but merely from the knowledge of the transverse section 
used to determine it. As we will see in our examples, this simplifies calculations considerably. 



In what follows we will assume that our manifold is homogeneous M = G/H, with H a closed 
subgroup. We denote by o £ G/H the distinguished class of H. We will finally show that given 
K s , and assuming p s ■ o = x s , for the moving frames generating K s , the A^-gon is completely 
determined up to the action of the group. 

Proposition 3.13. Let (x r ) and (x r ) be two twisted N-gons with left moving frames p and p 
such that p s ■ o = x s , p s ■ o = x s and p~ 1 p s+ \ = p~ l p s+ \ = K s . Then, there exists g E H such 
that x r = gx r for all r. 
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Proof. Notice that from p s+ \ = p s K s and p s +\ = p s K s it suffices to show that there exists 
h € H such that po = hp$. Indeed, if this is true we will have p s = hp s for all s, and from here 
p s ■ o = x s = hp s ■ o = h ■ x s . 

On the other hand, there clearly exists g £ H such that po = 5Pcb as it suffices to choose 
g = PoPq 1 , and g G H since both po and po leave o invariant. The proposition follows. □ 

Notice that the choice of o in this proposition is to some extent arbitrary (although it is geomet- 
rically the most convenient base point), and the assumption M = G/H is indeed not needed. 
One could have chosen any fixed point p G M, define p s • p = x s , p s ■ p = x s , and obtain g in the 
isotropy subgroup of p. One can even choose different base points for p and p in the same orbit 
and obtain a general element g £ G. 

4 Lie symmetric evolutions, maps and their invariant izat ions 

In this section we will show how to write any invariant time evolution of twisted A-gons, as 
well as any invariant map from the space of twisted A-gons to itself, in terms of the invariants 
and the moving frame, in a straight forward and explicit fashion. For simplicity we continue to 
assume that our manifold M is G/H and the normalisation equations include p s ■ o = x s (for a 
left frame) for all s. The analogous theorems in the continuous case were published in |34j . 

Definition 4.1. An evolution equation is said to have a Lie group symmetry if the Lie group 
action takes solutions to solutions. A recurrence map is said to have a Lie group symmetry if it 
is equivariant with respect to the action of the group. 

Such equations and maps are usually called invariant, and we will do so here. However, it is 
important to note that this does not mean that the equations are comprised of invariants of the 
action. 

Definition 4.2. We say the evolution 

(Xs)t = fs((Xr)) (14) 

is an invariant time evolution of the twisted A-gons under the action of the group G if its Lie 
symmetry group is G, that is, if (x r ) is a solution, so is (g • x r ) for any g G G. 

Denote by$ s : G/H — > G/H the map defined by the action of g G G on G/H, that is 
$ g (x) = g ■ x. Further, denote by T<3? 9 (z) the tangent map of <fr g at z £ G/H. The map T<fr g (z) 
is written as T<1> 9 if the base point z is clear from the contexlQ In local coordinates T$ ff is given 
by the Jacobian of <& g . We have by the chain rule both that (g ■ x s )t = TQ g (x s ) ((x s )t) and since 
$> g h = & g o we have that T^^^ = T<& g o T$^. A recent textbook reference on group actions 
and the associated multi- variable calculus, with and without co-ordinates, is [29] . 

Assume we have a discrete left moving frame p in the neighbourhood of a generic twisted A-gon, 
with a left action, so that p s {{g ■ x r )) = gp s ((x r )). We can now state the following Theorem. 

4 Given a manifold M and a map F : M — > M, one standard coordinate free definition of the tangent map 
TF(z) : T Z M -> T F(z) A/ is as follows. If a path j(t) C M satisfies 7(0) = z and the vector v G T Z M = 
d/dt| t=o7 (t), then TF(*)(v) = d/dt\ t=Q F ( 7 (i)). 
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Theorem 4.3. Any invariant evolution of the form \14y can be written as 

(x s ) t = Td> p »(v s ) (15) 

where o = [H], and v s ((x r )) G T Q M is a vector with invariant components, that is, v s ((g-x r )) = 
v s ((x r )) for any g £ G and for all s. 

Proof. We need to prove that 

T^^o)- 1 /^)) 

has invariant components under the action of G, and call it v s ; recall that the normalisation 



equations include p s • o = x s for all s, so that v s E T Q M. Indeed, since (14) is an invariant 
evolution, 

(g ■ x s ) t = T<5> g (x s ) ■ (x s ) t = T<S>g(x s ) • f s ({x r )) = f s ((g ■ x r )). 
Also, differentiating <& gPa = ^g^p 3 at o and applying the chain rule yields 

T$ Pa(ig . Xr)) (o) = T$ gM(xr)) (o) = T$ g (x s ) o T* p . ((av)) (o). 

From here 

T^ s((9 . :Cr ))(o)T^ 1 ((;E? . )) (o)/ s ((x r )) = fs(( 9 ■ X r )), 

which implies that TQ ps (( Xr ^(o)~ 1 f s ((x r )) is invariant. □ 

We now consider invariant maps. 
Definition 4.4. We say the map 

F(x s ) = (F s ((x r ))) , (16) 

is an invariant map of the twisted iV-gons under the action of the group G if its Lie symmetry 
group is G, that is, F((g • x s )) = gF((x s )) for any g £ G. 

We have the following theorem. We continue to assume we have a left moving frame and a left 
action. 

Theorem 4.5. If F is an invariant map of the form then 

F s ((x r )) = p s {{x r )) ■ z s ((x r )) (17) 
where z s {{x r )) G G/H is an invariant element, that is, z s ((g ■ x r )) = z s ((x r )) for any g G G. 

Proof. It suffices to show that pj 1 • F s ((x r )) is invariant and to call it z s . For a left moving 
frame and a left action we have 

Ps{(g ■ x r )) = gp s ((x r )) 
and by invariance of the map we have 

F((g ■ x s )) = (F s ((g ■ x r ))) = gF((x s )) = g (F s ((x r ))) , 

so that 

p-\{g ■ x r ))F s {(g ■ x r )) = /^((x^F^)). 

□ 
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If either a map or a time evolution is invariant, then there is a corresponding induced map or 
evolution on the invariants themselves. The reduction process is at times very involved and time 
consuming. Here we will describe a simple and straightforward way to find explicitly this so- 
called invariantization of the evolution. Further, in Section[5]we detail how, for a proper choice 
of evolutions, the resulting invariantizations are integrable, in the sense that they can be written 
in two different ways as a Hamiltonian system, using a Hamiltonian pair. We also illustrate, 



in Section 5.1.2, that some invariantized map also results in a discrete integrable (biPoisson) 
mapping. 

Assumption. From now on we will assume to have chosen q : M = H — > G, a section of the 
quotient G/H such that q{o) = e £ G, where e is the identity. 



Theorem 4.6. Assume we have an invariant evolution of the form (15) and let q be a section 
such that s(o) = e € G. Assume p r ■ o = x r and p r = ^(xr)pf 1 , where py £ H . Then 

(K s ) t = K S N S+1 - N S K S (18) 

where K s is the left Maurer-Cartan matrix and N s = p~ 1 (p s )t £ 0. Furthermore, if we split 
q = m © f), where q is the algebra of G, f) is the algebra of H and m is a linear complement 
that can be identified with the tangent to the image of the section and if N s = Ns + N™ splits 
accordingly, then 

N™ = T?(o)v s . (19) 

Proof. The first part of the proof is a straightforward computation 

(K s ) t = (pj 1 p s+ i) t = Ps l (p s+1 ) t - p~ l {p s )tp~ l ps+l = K S N S+1 - N S K S . 

Before proving the second part we notice that if ? is a section, then 

g^x) = q(g-x)h(x,g) (20) 

for some unique h £ H. In fact, one can take this relation as defining uniquely the action of the 
group G on a homogeneous space G/H in terms of the section. We will use this relation shortly. 

Since p s = s(x s )h((x r )), by the product and chain rules of differentiation, we have 

N s = p-\ Ps ) t = (TL prl ) (TR p? ) T,(x s )(x s ) t + (TL;i) Tpf ((x r ))((x r ) t ) (21) 

where L and R signify left and right multiplication and where TL and TR are their tangent 
maps. Very clearly the second term belongs to f) and so we will focus on the first term. If we 



differentiate ( 20 ) evaluated on x r we have 



(TL g )T<;(x r )(x r ) t = {TR h{Xr , g) ) T^ g (x r ))T^ a (x r )(x r ) t + {TL <g . x) ) Th(x r ,g)(c 



r)t 



for any g, where Th(x, g) is the derivative of h(x,g) as a function of x. If we now substitute 
g = p~ l and we use (15), ?(o) = e and T& p -i(x r ) = (T<J> Pi , (o)) _1 we have 



^ TL Pr 1 ) T< ;{x r ){x r )t = \ TR h ^ rtP -^J T?(o)v r + Th(x r ,p r ^(av)*- 
The last term is again an element of f) and we can ignore it. Finally, from ( |20[ ) we see that 

(pf 1 )~ 1 = p~ 1 q(x r ) = ■ X r )h(x r ,p~ l ) = sio^Xr^Pr 1 ) = h(x ri p~ l ). 
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Therefore 

Going back to the splitting of N r in (|21|) we see that 



N? = [TRpHj {TR ip?) -i) T,(o)v r = T,(o)v r 
as stated in the theorem. □ 



It is often the case that conditions (19) allow us to solve explicitly for the N r directly from 



equation (18), as we will see in the examples in the next section. Before going there we do a 



quick description of the invariantization of invariant maps. We will next prove the analogous 



result to Theorem 4.6 for invariant maps 



Theorem 4.7. Assume F is an invariant map given as in {lty . Extend this map naturally to 
functions of (x r ) using the relation F(l({x r ))) = £(F({x r ))). We will abuse notation and denote 
both maps with the same letter. Then the map induced on the invariants is given by 

F{K S ) = M~ l K s M s+1 (22) 

where M s = p~ 1 F(p s ). Furthermore, if ? is a section as before, and if p s = s(x s )pf for some 
pf E H, then M s = q{z s )M^ where Mf € H. 

Proof. First of all, notice that since by definition F(p s {(x r ))) = p s (F(x r )), F(p~ 1 )F(p s ) = 
F{p- l p s ) = F{I) = I and so F{p~ l ) = F{p s )~ l . From here 



F(K S ) = Fip-^Fips^) = F(p s )- 1 p s p- 1 p s+1 pjl 1 F(p s+1 ) = M~ x K s M t 



s+l- 



Also, assuming that p s = s(x s )p^ , a direct calculation using (20) shows 



M s = pJ 1 F(,(x s )p*) = pj\(F s ((x r )))F(p*) = pj\(p s .z s )F(p?) 
= p- l p s ^z s )h{z s , p s )F(pf) = <i{z s )h(z s , Ps )F(pf). 

It suffices to call = h(z s , p s )F(p^) G H to conclude the theorem. □ 



5 Completely integrable systems associated to discrete moving 
frames 

In this section we will describe invariantizations of general invariant evolutions for our previous 
two examples (centro-affine and projective) and we will associate completely integrable systems 
to each one of them by choosing the invariant elements v s defining the equation appropriately. 
In the centro-affine case we study the invariant maps and produce a well-known biPoisson map 
resulting from the invariantization of a particular invariant map among polygons obtained when 
choosing specific values for the invariant elements z s . Finally, we introduce a more complicated 
example, that of the homogeneous 2- sphere S 2 = SO(3)/SO(2). 
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5.1 Centro-affine case 



5.1.1 Invariant evolutions 



In the centro-affine case (see Example 3.9) the action is linear and so d& Ps (o)v s = p s v s . The 
general invariant evolution is given by 



\X S j X s -\-l | 



-x s+1 + v 2 s x s = p s 1 



where v].,v 2 are arbitrary functions of the invariants previously obtained in (11) (recall that p. 



was in this case a right frame and we need a left one here). If we want to find the evolution 
induced on kl and k 2 as in (11), then we will recall that the space H is the isotropy subgroup 



of e2, which is the subgroup of strictly lower triangular matrices. A section for the quotient is 
given by 

%- 1 







(23) 



if b 7^ (we work in a neighborhood of o = e2). Clearly p s = q(x r )p^ since ?(x r ) 1 p s € H. A 
complement m to f) is given by the upper triangular matrices. In that case 



and so 



N, 



i i 



-v 2 vl 



where a s is still to be determined. From equation (|18|) we have 

(*2)t 








-v 2 +1 vl +1 
a s+ i v 2 s+l 



k 2 s v 2 +i - k\a s+ i + v 2 k 2 



k i K S V S+1 



kW 



"1±1 

kl 



a s k 2 



V -^ + a s kl 




The entry (2, 2) of this system is given by = k\a s + (kl) ujn, which allows us to solve for 
the missing entry 



v l+i 



(k 



1\2- 



The other entries give us the evolution of the invariants. These are 



(k 2 s ) 



klv 2 +1 + k\v 2 



-k 2 v 2 



u s + l 

f- + k 2 v 2 + s x ? vl 



A 



Jk^ Us + 2 



,1 



(24) 



where the matrix difference operator is 

A = 

and where T is the shift operator Ta s = a s +i . 



k\(T+l) 
-k 2 (T-l) 



-k 2 T 



:T 2 



(25) 
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Let us define a diagonal matrix 



V 



(T-l)- l fc 




and compute the pseudo-difference operator 



Any I kKT+lXT-irX k\k 2 s 



(26) 



(27) 



fei / «- s ft s ' fc i y 
which is clearly anti-symmetric. We denote it by %{k\, k 2 ]. 

Theorem 5.1. The operator 1-L[k\,k 2 ], given by (27) is a Hamiltonian operator. It forms a 
Hamiltonian pair with Hamiltonian operator 



n [kl,k 2 ' 



k\ 
-kl 



Proof. Let us introduce the following Miura transformation 



k 1 ' 



kl 



Its Frechet derivative is 



D 



(Ps,q s ) 



s + l 



Is q- pi 



k 1 

K s+1 



VsTh 



1 j \ 1 j 

Under this transformation the operators Ho[kl, k 2 ] and H[kl, k 2 ] become 

HolPs,qs} = D {Ps , qs) Ho[k 1 s ,k 2 s ]Dl sjqs - I U />*<". -T) 



-(l - r-> 







and 



% s ,g s ]=I} (psigs) H[^,^ 2 ]£>* 



(Pb,9s) 



(28) 



Jfc 



o 1/ v - k Us 



T p s -PsT, 



PsT 1 p s - p s Tps p s (l - T)q s 
-q s (l- T' 1 )p s T~ x p s - p s T 

These two operators form a Hamiltonian pair for the well-known Toda lattice Q in Flaschka 
coordinates [H[39]. Indeed, we have 

and thus we proved the statement. □ 



becomes 



If we take ( v 2 , Vg) = (0, —kl_ 1 ) in (24), then the evolution of the invariants for kl and k 2 



4)t 



k 1 

K s-1 



k x k 2 



-.(k 2 ) 2 k 1 



k 1 



%5ki 



{k s )t I \ ti 1 ,,i s I V 2 

which is an integrable differential difference equation and can be transformed into Toda lattice 



(29) under the transformation (28). We summarize our results as a Theorem. 
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Theorem 5.2. The evolution of polygons in the centro-affine plane described by the equation 

k 1 

(X s )t = -, rX s+ l 



\X S , X s -\-l I 

induces a completely integrable system in its curvatures k\,k 2 equivalent to the Toda Lattice. 



The pseudo- difference operator V is really a formal operator as T — 1 is not invertible in the 
periodic case. It means that we apply the operator only to Hamiltonians whose gradients are in 
the image of T — 1, as is the case here. 

5.1.2 Invariant maps 

In the centro-affine case an invariant map is of the form 

7-1/ \ z s 2 —1 ( %\ 

F(x s ) = -- -x s +i + z s x s = p s | 

\%Sl %s + l | \ z s 



Although in this simple case we could directly find the transformation of the invariants, we will 
follow the process described in Section |4j Our example in Section 5.3 will provide a stronger 



case for the effectiveness of the method. Using section (23), the transformation under F of kl 
and kl will be built up using the matrix M s = — p s+ ip^^given by 

M _ ( fe 2 )" 1 z\\ ( 1 

where a s needs to be found. The equations relating K s and M s are given by 
k 2 -k ls 



F 



4 -4 \ ( kl -h\\ A^)" 1 ( i o 

The (2, 2) entry of this equation will allow us to solve for a s . It is given by 

(1 1 1 ,2/2 7,1 1 £,2\\ , Z s+1 _ n 
TJ Z s Z s+l "+~ z s\ z s+l K s z s+l K s) I 2TT ~~ 

K s J Z s K s 

The other entries of the system will solve for the transformation of kj and k 2 . They are given 
by 

^(^s) = jA Z s Z s+l + ( z s+l^s ~ z l+l^) z s 
K s 

T?( 1 2 \ ( f 7 2 1 2 7 1 \ 2 1 11 \ -^s Z< s 

* K^s) = I \k s z s+l ~ z s+l fi s) z s ~~ T\ Z s z s+l I a s+l 2 El 2 

V K s / Z s+l K s z s+l 

Now we take z 2 = ^, where c 7^ is constant and z\ satisfying the relation 

1 k 1 z 1 t k 2 

TT z s z s+i "I f" = z s+i (30) 
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Then the above maps become 



F(kl) 
F(k 2 s ) 



's+1 



'S+1 



K s+1 S 



+e 



Let a, 



K s + 1 



kg and f3 s 



Then using (31) and (32) we have 



F(a s 
F(b s 



F(kl) _ k] (3 S+1 _ f3 s+1 



F{k\ +1 ) (3 s kl +1 s P. 
k\ k s _i 



b s + c 



k\ + ifis k\fi s -i 



b s + c 



a s a s -i 



Ps Ps-1 



The constraint (30) on z\ becomes 



a,s-i , 1_ 



that is, 



Ps-lPs 



p s = l + cb s 



1 



CPs Ps 



2 a s-l 



P. 



8-1 



(31) 
(32) 



Thus we obtain the integrable discretization of the Toda lattice as the formulas (3.8.2) and 
(3.8.3) in [39] . Thus we obtain the following result: 

Theorem 5.3. The invariant map of polygons in the centro-affine plane described by the equa- 
tion 



F(x s ) 



I rXs+i H 

\^8 5 ^8+1 | ^ 



where c ^ is constant and z\ satisfying (30) induces a completely integrable map in its curva- 
tures k~,kg equivalent to the integrable discretization of the Toda Lattice. 



5.2 The projective case 

In the projective example the subgroup H is the isotropy subgroup of 0, that is, H is given by 
lower triangular matrices. Thus, a section can be chosen to be 



1 x 
1 



One can check directly that <;{x r ) l p r G H and also 

*io)v=(» i r 

A general invariant evolution in this case would be of the form 

(x s ) t = d$ Ps (0)v s = \v s 



(33) 



23 



where v s is a function of (k r ) and where a s is the (2,2) entry of p s . The evolution induced on 



k s is given by the equation ( 18 ) with 



Ps ~OL t 



We have 

'{k s )t 0\ _ ( k s l\ ( a s +i v s+ i \ (a s v s \ ( k s 1 



OJ \-l OJ \P S+1 -a s+1 J \p s -a s J V-l y 

The entries (2, 2) and (1, 2) of this system will allow us to solve for f3 r and a r . Indeed, the (2, 2) 
entry is given by 

= -P s ~ v s+ i 

and so 

Ps = -V s+ i. 

The (1, 2) entry is given by 

a s+ i + a s - fc s v s+ i = 0. 

This leads to 

a s = (T+ l) _1 A; s v s+ i, 

where T is the shift operator. Notice that T + 1 is an invertible operator. Assuming N is not 
even we can solve the equation (T + l)a s = w s to obtain 

\i=0 i=s ) 

Finally, the (1,1) entry gives the evolution for k s 

(k s ) t = v s - v s+2 + k s (a s+ i - a s ) 

= (r- l -r+k s (r-i)(r+i)- 1 k s )v s+1 . (34) 

Theorem 5.4. The anti- symmetric operators 

H 1 [k s ]=T-T- 1 

and 

n 2 [k s ] = k s (T-l)(T+l)- 1 k s 

form a Hamiltonian pair. 

Proof. Let us introduce the following Miura transformation 

u s = 7T • ( 35 ) 

Ks 



We compute its Frechet derivative D/ is = —\ju 2 s . Under the transformation (35) the operators 
Hi[k s ] and H2[k s ] become 

% x [us] = D^-H^D*- 1 = u 2 s (T-T- 1 ) u 2 s - 

H 2 [u s ] = D k ^n 2 [k s ]Dl- 1 = Us (T-l)(T+ l)~ l u s . 

These two operators form a Hamiltonian pair for the well-known modified Volterra lattice Q 



(u s ) t = u 2 s (u s+1 - u s _i) = tii[u s ]5 Us lnu s = U2[us\8 Ua {usU s -i)- (36) 
This implies that Hi[k s ] and %2 [k s ] form a Hamiltonian pair and the statement is proved. □ 
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Take v s +i = 5k s lnk s in (34). The evolution for k s becomes 

1 



1 



sit 



ks 



k 



s+1 



which is an integrable difference equation. Under the transformation (35 ), it leads to the modified 



Volterra lattice ( 36 ) . 



At this moment it seems natural to wonder if one could also obtain this modified Volterra 



evolution in the centra- affine case (5.1), via a reduction to the space k\ = 1 for all s. And 
indeed such is the case. 

Assume k\ = 1 and choose those evolutions that leave k\ invariant. That is, assume (k\)t = 0. 

,2 



Then, from (24) we have 



v s+i + v « ~ k^ yl 



s u s+l 







and from here 



Substituting in ( 24 ) we get 



(k 2 s )t 



Hence, 



-k 2 s T{T+ l)- l k 2 s Tv] - v\ + k 2 s {T +l)- l k 2 s Tvl + 7"V 

= [T-T- 1 - k 2 (T- 1)(T + l)- 1 ^ 2 ] Tvl 

which is identical to the projective one with identifications k s = k 2 and v s 
choosing Tv\ = — 5^.2 In | A; 2 1 will result in a modified Volterra equation for the centro-affine 
invariant A; 2 , as far as the initial polygon satisfy k\ = 1, that is, as far as the area of the 
parallelogram formed by x s and in the plane is equal 1 for all s and they are properly 
oriented. Other constant values can be chosen with minimal changes. Notice that this evolution 
is well-defined in the periodic case. 



5.3 The homogenous sphere S 2 SO(3)/SO(2) 

In this section we consider invariant evolutions on the homogeneous sphere, S 2 = SO(3)/SO(2). 
We first consider a local section using which we can describe our calculations of a discrete frame, 
the associated discrete Maurer-Cartan matrices, and invariantizations. We then show that, for 
a certain choice of polygon evolution the resulting curvature flow is integrable, of Volterra type. 

If G = SO (3) and H = SO(2), we consider the following splitting of the Lie algebra into 
subspaces so (3) = m © f) with 



y\ (A 

o em 



t n J e m n jGt, (37) 



where y G M 2 and A £ so(2). Associated to this splitting we have a local factorization in the 
group into factors belonging to H = SO(2) and exp(m). This factorization is given by 

O 0\ f I + cosy yy T sin y y 



g = g(Q, y) = " " ' ^ « (38) 

yv ,yj y Q 1 j y - smy yT cos(\\y\\)J v 1 

, cos(||y||) — 1 . sin(||y||) 9 T . . .in. 

where cos„ = , sm„ = — r—r — and \\y \\ = y y. ihe factorization exists locally m a 

\\y\\ Nl 

neighborhood of the identity. 
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Let ? : S 2 — > S0(3) be the section denned by the exponential, that is 

/ \ _ ( x\ _ (l + cos x xx T svn x x 
?(Xj " eXP V-x T oJ"V-sin x x T cos(||x||; 

where we are using local coordinates around the south pole (whose coordinates are zero). One 
clearly has that tfc(o) : T Q S 2 — > m is an isomorphism given by 

*(o)«=(_Jt o 

The action of SO(3) on the sphere associated to this section, let's denote it by g-x, is determined 
by the relation 

gq(x) = s(g • x)h 



for some h E SO(2) which is also determined by this relation. Let g be as in (38). Straightforward 
calculations show that, if rj = g ■ x, then 

sin^ r\ = sin x Qx + (cos^ sin x y T x + sin y cos(||x||)) Qy (39) 

and 

cos (1 77!) = cos (1 7/1) cos(||x||) - sm x sm y y T x. (40) 

The last equation can be expressed in terms of the cosine of a certain angle using the spherical 
law of cosines, as we will see later. To have a better idea of what the coordinates given by this 
section are, recall the standard geometric identification of SO(3)/SO(2) with the sphere: 

Given an element in SO (3) we can identify the last column with a point on the sphere and the 
first two columns as vectors tangent to the sphere at the point. The element of SO (2) acts on 
the two tangent vectors. With this identification, our coordinates result on 



P 



smx H 
cos (Ixl) 



being the point on the sphere. If we consider 9 and 4> to be the standard spherical angles, then 

||x|| = 4> an d |f| = ^g-^^- Therefore, the coordinates x describe the projection of p on the 

xy plane, multiplied by the angle <j>. See the picture below. This might seem as a cumbersome 
choice, but the advantages in calculations will be worthy, plus our Serret-Frenet equations will 
look very familiar to the reader. 



5.3.1 Moving frame and invariants 

Next we will determine a moving frame using the normalization constants c r T = and c r r+1 = a r e%, 
where a r is an invariant still to be identified. Assume p r = g(6 r ,7/ r ) as in (38). The equation 

determines the choice y r = —x r , and p r ■ x r+ \ is determined by the equation p r ■ x r +i = a r e\ 
implicitly written as 

sin ar£l a r ei = <d r (sin Xr+1 x r+ i + (cos Xr sm Xr+1 x^x r+1 - sin Xr cos(||x. r+ i|))x r ) . (41) 

The invariant a r is determined by the condition O r G SO(2), while this equation determines O r . 
If we impose the further condition a r > for all r, then O r is uniquely determined (it might 
belong to different connected components of SO(2) for different r's). 
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X r +1 

Figure 4: Local coordinates on the sphere 



Theorem 5.5. The right Maurer-Cartan matrices associated to p r as above are given by 

I k r 0\ /0 -a r \ 
JT r = exp \-k r exp . (42) 





We caZZ a r the discrete spherical arc-length invariant, and k r the discrete spherical curvature of 
the polygon. 



Proof. Let us call K r = g(T, k) as in (38). To prove the theorem we can use the recursion 
equations (13) for K r to determine it. For a right matrix the equations are given by K r ■ c£ , 1 = 
c r+i) which, when substituted in (39), and after minor simplifications, becomes 

= sin ar a r e\ + (cos K sin ar a r n T ei + sin K cos a r ) k. 

Let us re-write this as = sin ar a r e\ + X r n so that k = — sul y e\. We can in fact use the 
expression for X r to solve for it. Indeed, again after minor simplifications we obtain 



X r = cos K sin ar a r n T ei + sin K cosa r = X r I 1 



cos 



sin a r 
X r 



cos a r 

H sin 

sin a r 



sin a r 
X r 



This results in the relation 
' sin a r 



sin 



X r 



I sin a, 



cos a r 



cos 



which results in X r = sin a and k 



a r ei. 



sin a r 



sin 



sin a r 

X r 



This calculation determines the m component of K r as in the statement of the theorem. Notice 
that the I) component T depends on only one parameter and therefore we can write it as in the 
statement, even if we do not give its explicit formula. □ 



Both invariants a r and k r have a very simple geometric description, as shown in our next theorem. 

Theorem 5.6. The spherical arc-length invariant a r is the length of the arc joining x r to x r +i- 
Let f3 r ^ s be the counterclockwise angle formed by the arc x r x s and the arc x s N , where N is the 
north pole. Then, 

K-l = VT — (P r ,r+1 + Pr+l,r+2) 
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Figure 5: Description of the discrete curvature. Notice that when /3 r+ i ir+ 2 + Pr,r+i = ^ the 
curvature is zero. 



Proof. Calculating the length of the vector in equation (41) we see that 



sin 2 (a r ) = sin ||x r +i|| +((cos||x r || — 1) sin |av+i| cos 9 r ,r+i — sin ||x r || cos ||x r +i||) 

where # r ,r+i = r — # r +i and 9 r is the polar angle of x r . After some simplifications this equality 
becomes 

sin 2 (a r ) = 1 — (cos ||x r || cos + sin \x T \ sin ||x r +i|| cos r ,r+i) 2 ■ 

Using the spherical law of cosines (see figure 3) we obtain that 

cos \\x r \\ cos ||x r +i|| + sin ||x r || sin |x r +i|| cos(6> r ,r+i) = cos a rt r+i, 

where a r , r +i is the angle between p r and p r +i, the points on the sphere corresponding to x r and 
x r+ i. Using this fact, the above becomes 

— cos a T: r + i + 1 = sin a r ^ r+ i = sin a r . 

Therefore, a r = \a r>r +i\ since it is a positive number, as stated. 

Law of cosines: 

cos c = cos b cos a + sin b sin a cos C 

Law of sines: 
sin a sin b sin c 




sin A sin B sin C 



Figure 6: Spherical trigonometric law of cosines and sines. The angles A, B and C are the 
corner angles while a, b, c are the angles corresponding to the indicated arc as measured from 
the center of the sphere. 

A similar calculation, although a more involved one, describes k r . Using the recursion relation 

A f L r+1 — c r+l' 
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the fact that c^i = 0, the expression of K r and the formula (39) for the action, we have that 



sin r+2 eft? = sin a r 

S-+1 '' l 



- cos 
sin k r 



Therefore, tt — k r is the polar angle of c^t\ (while fc^il = a r ). Now, we do know that p r+ 2 
x r+ \ = cjlijj where p is the right moving frame; therefore, we can again use (39) to find 

„r+2 



sin r+2 c \\ — Q r+ 2V r+ i r+ 2 

S-+1 ~ 



(43) 



where 

Vij = sin Xi Xi + (cos Xj sin Xi xj ■ x« - sm Xj cos ||x;||) Xj. 

We clearly see that the polar angle of c£i 2 is a combination of the polar angle of Q r +2 an d 
the one of V r+ i^ r+ 2- Now, we know that r is determined by the equation (41). That is, ® r +2 
rotates y r _|_3 )r+ 2 to the x-axis 

9 r+2 (sin :Cr+ 3 x r+3 + (coSa; r+2 siHa. r+3 x r+2 ■ x r+3 - sin Xr+2 cos ||x r+3 |) x r+2 ) 
= @r+2Vr+3,r+2 = sin ar+2ei a r+ ie\. 
If we denote by P(v) the polar angle of the vector v, from here we get that 7r — k r = P(V r +x,r+2) — 
P(Vr+3,r+2)- 

Next, we can write x r in terms of the polar coordinates of p r , that is, 



sin \\x r \\ cos v r 
sin I x r I sin 6 r 



and substitute in the expression for Vij. After some minor trigonometric manipulations we 
obtain 



Vi 



cos 
sin 



sin ( 



cos ( 



Form here P(V r +i 



r+2) 



' cos \\xj\\ sin \\xi\\ cos(6*jj) — sin \\xj\\ cos \\xi 
sin ||xj|| sin(6ij) 

P(V r +3,r+2) = P{v r +i,r+2) - P(v r+3)r+2 ), where 

^cos \\xj\\ sin ||xj|| cos(^jj) — sin \\xj\\ cos \xi 



sm xj sin 



"1,3; 



(44) 



Notice that, as before, after basic trigonometric manipulations 



1 



cos 



(||xi||) cos(|xj|) +sin(||x i ||)sin(||x : ,-||)cos(6'jj)) 2 = 1 - cos 2 a itj 



sm a 



where we have used the spherical law of cosines once more. Notice that < aij < tt and so 
from (44) the sine of the polar angle of Vij is given by 

sin ||xj|| sin(#jj) 



Let us assume all the angles involved are in the first two quadrants. Then, the spherical law of 
sines tells us that 

sina r +i ir +2 sin||x r || 
sin# r +i,r+2 sin/3 r+ i ;7 . + 2 
and so the sine of the polar angle of t; r+ i ir+ 2 is given by 



sm ||x r+ i|| sm^ev+i^-^j 



sin/3, 



r+l,r+2- 



sm a r +i, r +2 

Likewise, the polar angle of f r +3, r +2 is /3 r+ 3 ir+ 2 = — (3 r+ 2, r +3- Therefore, 

TT — k r = (3 r+ i ir+ 2 + /3 r+ 2, r +3 

as stated. 



□ 
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5.3.2 Invariantization of invariant evolutions and its associated integrable system 



We will next describe the evolution induced on K r by an invariant evolution of polygons on the 



sphere. According to (15) and using our definition of moving frame, the most general form for 



an invariant evolution of polygons on the sphere is given by 

T \ T 

II'''' 

fr l \ 

for some invariant vector r s = I | I depending on (k r ) and (a r ). 

Theorem 5.7. Assume the twisted polygons (x r ) are solutions of an invariant evolution of the 
form above. Then the invariants a s and k s evolve following the equations 

(a s ) t = -t\ + r] +1 cos k s - r 2 s+l sin k s (45) 

ow, = r , i _ r , «^ +r; ^ +r , ^ (46) 

sma s smo s +i sina s +i sma s +i 

l cosa s 2 cosa s 

- r s+1 - smfc s -r s+1 - cosk s (47) 

sma s sinoj 

Proof. Assume the polygons are evolving following this evolution; then, our left Maurer-Cartan 
matrices K s = K~ l will follow the evolution 

K- l (K s ) t = N s+1 - K^NsKs (48) 



where 






n s 




-n s 







-r 1 

' s 


-r 2 

' s 






Next, we remark that I T K I commutes with its derivative since n is a multiple of e±. This 

. K U 



means 



o -«\\ / ((Is) M fo -k 
( ' xp U T o JJ,= I ° J ° ) exp U o 



Using the expression of K r in the previous theorem, one has that 

-(k s ) t (a s ) t cos k s 

K- l {K s ) t = \ (k s )t -(a s ) t smk 

-(a s ) t cos k s (a s ) t sink s 



Substituting this into ( 48 ) and multiplying out matrices we get 

(k s ) t = -n s+ i + n s cos a s + r 2 s sin a s 

and 



cosk s sh\k s \ ({a s )t \ _ ( cosk s sink s \ ( —r s \ t , s+l 



sm.k s cosk s J V y V — sin& s cosk s J V— r 2 cosa s + n s sina Sy / ^ \^r 2 +1/ 



„l \ / r i 
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From here we get 

2 cosa s 1 , 1 .7.2 ; \ 

n s = r s ; (r s+1 sin k s + r s+1 cos k s J 

sma s sma s ' 

and 

(a s ) t = -r\ - r 2 s+l sin k s + r^ +1 cos k s . 
Substituting the vale of n s in the evolution of k s completes the proof. □ 

In particular, if we ask that the evolution preserves the spherical arc-length (i.e. (a s )t = 0) and 
we choose the value tana s = 1 as constant value, then choosing r\ = 1 for all s produces the 
curvature evolution 

= ^ ( 1 -<***«+! _ l-COB^-A 

y sin/c s+ i sinfc s _i y 

Let u s = 1 ~^ fca . The above equation becomes 

/2 

(4 = y ( 49 ) 



which is a special case of equation (VI) in the list of integrable Volterra-type equations [42] . 
It is a bi-Hamiltonian equation [33], where the second Hamiltonian operator was not explicitly 



given. Here we give its compatible Hamiltonian and symplectic operators [5] . Equation ( 49 ) can 
be written as 

(u s ) t = —U5 Us ln(l + u 2 ) and Z{u s ) t = — S Ua ((l + u 2 )u s+1 u s -! + --u^Li) 

where 

n = (i + u 2 s ) (r-r- 1 ) (1 + ^) 

is Hamiltonian and 

1 = T-T- 1 + t^(T - l)" 1 ^! + u a -!) + K+i + n s -i)T(T - I)" 1 - 2 "~ 



1 + u% 1+U 



s 



is a symplectic operator. Furthermore, these two operators satisfy 

ut = k 2 = f (i + u 2 s )t + 2u s u s+l + (i + u^r- 1 + kmt - ir 1 :^^ 

V 1 + ut 



where K is a Nijenhuis recursion operator of (49). Its recursion operator K can not be written as 



the product of weakly nonlocal Hamiltonian and symplectic operators. A similar example was 
presented in |41j . 



Notice that, as before, T — 1 is not invertible in the periodic case. That means we will need 
to work on infinite gons, or assume we are working with Hamiltonians whose gradient is in the 
image of the operator, as it is the case with this example. 
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6 Conclusion and Future Work 



In this paper, we have developed a notion of a discrete moving frame and shown it has com- 
putational advantages over a single frame for the invariantization of discrete evolution flows 
and mappings. Further, we have shown in our examples that the use of discrete moving frames 
greatly aids the identification of discrete integrable systems and biPoisson maps that can be 
obtained as invariantizations of invariant flows of gons. 

Investigations are under way on how the Hamiltonian structures might appear in the general 
case, and in particular how they appear in the projective plane. Further work also remains to 
illustrate and detail how our constructions of invariant mappings and their invariantizations lead 



in general to discrete integrable mappings, as seen in Section 5.1.2 This work is different in 
nature to the previous one and it will be greatly aided by understanding the differential-difference 
problem. 

Hamiltonian structures in the projective plane could be relevant to the results in [38]. There 
the authors proved that the so-called pentagram map is completely integrable, but they did not 
provide a biHamiltonian structure for its invariantization. A structure obtained in this setting 
would be a natural candidate to prove that the pentagram map is biPoisson. 

The study of discrete moving frames for more general applications, and insight into how our 
theorems in Section 3 here may generalize, remains to be achieved. One such application is to 
the discrete Calculus of Variations, as briefly indicated in Example |2.9| We believe however, 
that eventually the applications for a more general theory of discrete moving frames will extend 
to invariant numerical schemes and computer graphics, amongst others. 
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